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SENSITIVITY  ANALYSIS  OF  A  NONLINEAR 
STRUCTURAL  DESIGN  PROBLEM 

by 

Anthony  V.  Fiacco 
Abolfazl  Ghaemi 

1.  INTRODUCTION 

The  purpose  of  the  current  paper  is  to  conduct  a  first  order 
sensitivity  analysis  of  a  nonlinear  structural  design  pro^xem  discussed 
in  Chapter  6  of  the  book  by  Bracken  and  McCormick  [8]. 

By  sensitivity  analysis  is  meant  an  analysis  of  the  effect  on  the 
optimal  objective  function  value  and  on  an  optimal  solution  point  of  small 
perturbations  in  the  <model  parameters.  The  jjmportance  of  such  an  analysis 
in  real  world  optimization  problems  cannot  be  overstated.  It  provides  the 
model  maker  and  user  with  invaluable  information  regarding  the  functional 
relationship  between  a  solution  and  the  design  parameters.  This  has  many 
potential  applications.  For  example,  identification  of  those  parameters 
having  the  most  significant  impact  on  the  optimal  solution  can  provide  a 
basis  for  developing  educated  guidelines  for  taking  appropriate  and  efficient 
action  toward  effecting  parameter  changes  that  will  give  an  optimal  marginal 
improvement  of  system  design  or  performance. 

A  theoretical  basis  for  sensitivity  analysis  for  nonlinear  programming 
was  given  by  Fiacco  and  McCormick  [13]  and  generalized  and  extended  by 
Fiacco  [10].  Based  on  the  approach  given  in  [10],  an  algorithmic  procedure 
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was  proposed  by  Fiacco  and  implemented  by  Causey  |9).  Later  on  it  was 
refined  and  recoded  by  Mylander  [14],  using  the  several  subroutines  of  the. 
SlMT-version  4  computer  code  by  Mylander,  Holmes,  and  McCormick  [15], 
which  implements  the  Sequential  Unconstrained  Minimization  Technique  for 
nonlinear  programming  using  the  logarithmic-quadratic-loss  penalty  function. 
This  will  be  described  in  Section  2.  Mylander's  sensitivity  subroutines 
were  integrated  with  the  SUMT-version  4  computer  program  by  Armacost  and 
Mylander  [7|.  This  routine  was  further  revised  and  expanded  by  Armacost  [1]. 

The  latter  version  of  the  routine,  now  called  "SENSUMT,"  is  compiled 
in  The  George  Washington  University  Computer  Center  and  constitutes  the  main 
element  on  which  the  present  sensitivity  analysis  is  based.  The  study  is 
conducted  in  three  phases:  Phase  1  deals  with  the  solution  of  the  problem. 
Phase  2  deals  with  the  solution  and  sensitivity  analysis  study  of  the  problem 
when  each  of  the  right  hand  sides  of  the  constraints  is  perturbed.  Several 
problems  are  analyzed,  each  associated  with  a  significantly  different  right 
hand  side  initial  value.  Phase  3  deals  with  a  sensitivity  analysis  study 
with  respect  to  each  problem  design  parameter,  again  for  several  problems, 
each  associated  with  a  significantly  different  set  of  problem  design  parameter 
initial  values. 

Estimates  are  given  of  the  optimal  solution  of  this  family  of  problems, 
using  the  optimal  solution  and  first  order  sensitivity  information  of  the 
unperturbed  problem.  Furthermore,  several  stability  and  convergence  character¬ 
istics  of  the  solution  points  and  their  partial  derivatives  with  respect  to  the 
problem  parameters  are  computationally  verified. 

As  an  example  of  the  numerous  inferences  that  can  be  made,  the  analysis 
reveals  that  the  weight  of  the  optimally  designed  structure  is  most  sensitive 
to  the  minimum  allowable  structural  thickness  and,  next,  to  a  corrosion 
allowance,  two  of  the  many  design  parameters  of  the  problem. 

Similar  studies  have  been  conducted  by  Armacost  and  Fiacco  on  a  variety 
of  problems,  including  a  cattle  feed  problem  [2]  and  a  multi-item  inventory 
problem  [5];  and  by  the  authors  on  stream-water  pollution-abatement  model 
involving  numerous  parameters  [11],  [12]. 

-  2  - 
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2.  BASIC  SENSITIVITY  RESULTS 

A  comprehensive  review  of  the  sensitivity  results  obtained  for  nonlinear 
programming  problems  is  outside  the  scope  of  the  present  paper.  In  the  follow¬ 
ing  we  review  only  the  theory  supporting  the  validity  of  the  computational 
results  obtained.  The  basic  result  was  given  in  Fiacco  and  McCormick  [13, 
Theorem  6]  for  a  particular  class  of  perturbations.  Fiacco  generalized  the 
theory  and  also  established  a  theoretical  basis  for  utilizing  a  penalty  func¬ 
tion  method  to  estimate  sensitivity  information  associated  with  a  local  solu¬ 
tion  and  its  associated  optimal  Lagrange  multipliers,  for  a  large  class  of 
nonlinear  programming  problems,  with  respect  to  a  general  parametric  variation 
in  the  problem  functions  [10]. 

Armacost  and  Fiacco  subsequently  computationally  implemented  this 
approach  to  demonstrate  practical  applicability  [2]  and  applied  this  theory 
to  obtain  the  first  and  second  order  sensitivity  results  for  the  optimal  value 
function,  deriving  formulas  for  the  efficient  calculation  of  the  sensitivities 
of  this  function,  as  well  as  the  sensitivities  of  the  local  solution  point  and 
its  associated  optimal  Lagrange  multipliers  [3],  [4]. 

The  parametric  mathematical  programming  problem  considered  by  Fiacco 
is  of  the  following  general  form: 

minimize  f(x,e) 

x  e  E 

subject  to  g^x.e)  ^0  ,  i  =  l,...,m  ,  P(e) 

hj(x,E)  "  0  *  3  =  1 . .  » 

where  x  is  the  usual  vector  of  variables  and  E  is  a  k-component  vector  of 
numbers  called  "parameters."  It  is  desired  to  analyze  the  behavior  of  a  solu¬ 
tion  vector  x(e)  and  the  optimal  solution  value  f*(e)  £  f[x(e),e]  near 
some  given  value  of  c  .  Without  loss  of  generality*  assume  that  the  parameter 
vector  of  interest  is  e  *  0  . 

The  Lagrangian  for  Problem  P(c)  is  defined  as 

m  p 

L(x,u,w,e)  £  f(x,e)  -  E  u  g. (x,e)  +  E  w  h  (x,e)  . 

i-1  11  j-1  3  3 

-  3  - 
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The  sensitivity  results  are  based  on  the  following  four  assumptions: 

A1  -  The  functions  defining  Problem  P(t)  are  twice 

continuously  differentiable  in  (x,e)  in  a  neigh¬ 
borhood  of  (x*,0). 

A2  -  The  second  order  sufficient  conditions  for  a  local 

minimum  of  Problem  P(0)  hold  at  x*  with  associated 
Lagrange  multipliers  u*  and  w*  . 

A3  -  The  gradients  Vxg^(x*,0)  ,  for  all  i  such  that 

g  (x*,0)  =  0  ,  and  V  h  (x*,0),  j  =  l,...,p  are 
i  x  1 

linearly  independent. 

A4  -  Strict  complementary  slackness  holds  at  x*  when 

C  =  0  (i.e.,  u*  >  0  for  all  i  such  that  g^xfc.O)  =  0)  . 

Under  the  above  assumptions,  Fiacco  [10]  established  the  following  generalization 
of  Theorem  6  in  [131. 

Lemma  2.1  (Local  characterization  of  a  Kuhn-Tucker  triple.)  If 
assumption  Al,  A2,  A3  and  A4  hold  for  Problem  P(e)  at  (x*,0),  then 

(a)  x*  is  a  local  isolated  minimizing  point  of  Problem  P(0) 
and  the  associated  Lagrange  multipliers  u*  and  w*  are 
unique; 

(b)  for  r.  in  a  neighborhood  of  0  ,  there  exists  a  unique 
once  continuously  differentiable  vector  function 

T 

y(e)  =  (x(e)»u(e),w(e))  satisfying  the  second  order 
sufficient  conditions  for  a  local  minimum  of  Problem  P(t.) 

T 

such  that  y(0)  =  (x*,u*,w*)  =  y*  and  hence,  x(e)  is 
a  locally  unique,  local  minimum  of  Problem  P(e)  with  asso¬ 
ciated  unique  Lagrange  multiplers  u  (  e  )  and  w(c)  ;  and 

(c)  for  c  near  0  ,  the  set  of  binding  inequalities  is 
unchanged,  strict  complementary  slackness  continues  to 
hold,  and  the  binding  constraint  gradients  are  linearly 
independent  at  x(e)  . 
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(d)  (Fiacco  and  Armacost  [3]),  for  £  near  0  ,  the  gradient 

of  the  optimal  value  function  is 

V£f*(e)  =  V£L(u,c)  at  y  =  y(e)  , 

(e)  which  also  means  that,  for  £  near  0  ,  the  Hessian  of 

the  optimal  value  function  is 

V^f*(e)  =  V*L(y(£),e)  . 

The  above  results  provide  a  characterization  of  a  local  solution  of 
Problem  P(e)  and  its  associated  optimal  Lagrange  multipliers  near  e  =  0  . 
They  show  that  the  Kuhn-Tucker  triple  y(e)  is  unique  and  well  behaved, 
under  the  given  conditions.  Since  y(£)  is  once  differentiable,  the  par¬ 
tial  derivatives  of  the  components  of  y(E)  are  well  defined.  This  fact 
and  Assumption  A1  also  mean  that  the  functions  defining  Problem  P(e)  are 
once  continuously  differentiable  functions  of  £  along  the  "solution  tra¬ 
jectory"  x(e)  near  £  *  0  ,  and  the  Lagrangian  is  a  once  continuously 
differentiable  function  of  e  along  the  "Kuhn-Tucker  point  trajectory." 

The  above  results  constitute  the  structure  for  numerous  developments  and 
extensions,  many  of  which  have  been  established  by  Fiacco  [10]  and  Armacost 
and  Fiacco  [2-6]. 

The  realization  of  this  theorem  for  the  parametric  right  hand  side 
problem  of  special  interest  in  the  present  study  is  treated  in  detail  by 
Armacost  and  Fiacco  [' ].  The  parametric  right  hand  side  problem  is  the 
following  important  realization  of  P(£): 

minimize  f (x) 

subject  to  g^x)  _>  e  i=l,...,m  R(t) 

hj(x)  =  ej+m  i  =  1 . .  ’ 

The  Lagrangian  for  R(e)  is 

m  p 

L(x,u,w)  =  f(x)  -  Z  u  (g  (x)-e  )  +  Z  w. (h, (x)-e.  )  . 

i-1  11  1  j»i  J  J  J+m 

As  evident  from  the  Lagrangian,  the  results  (d)  and  (e)  of  Lennna  2.1  for 
problem  R(e)  simplify  respectively  to 
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(d*)  Vf  f*(>) 

(e')  V*  f*(r.) 

Fiacco  [10]  has  shown  that  the  class  of  algorithms  based  on  twice 
continuously  differentiable  penalty  functions  -(specifically,  using  the 
logarithmic-quadratic  loss  penalty  function)  can  be  used  to  estimate  y(>  ) 
and  its  derivatives  in  a  neighborhood  of  t:  =  0  ,  for  the  general  problem 
P(t)  .  Minimization  of  the  penalty  function  with  penalty  parameter  r 
yields  a  solution  of  a  perturbation  of  the  ICuhn-Tucker  system  in  a  neighbor¬ 
hood  of  (t  ,r)  =  (0,0)  .  Armacost  and  Fiacco  [3]  define  an  optimal  value 
penalty  function  and  obtain  first-  and  second-order  sensitivity  estimates 
which  converge  to  the  corresponding  sensitivities  for  the  optimal  value 
function  for  Problem  P(t). 

The  logarithmic-quadratic  penalty  function  is 

m  P 

W(x,f  ,r)  f(x,r.)  -  r  I  £n  g.(x,r)  +  (l/2r)  7  h7(x,f)  . 

i-1  j=l  1 

l.orana  2 . 2  (Fiacco  flO,  Theorem  3.11).  If  the  assumptions  Aj  ~ 
hold,  then  in  a  neighborhood  of  (e,r)  =  (0,0)  there  exists  a  unique  once 

T 

continuously  differentiable  vector  function  y(f  ,r)  =  [x(t,r)  ,u(r.,r),w(f  ,r)  ] 
satisfying 

VxL(x,u,w,c)  =  0 

u^g^(x,e)  =  r  ,  i=l,...,m  ,  (2.1) 

hjCx.e)  -  vy  ,  j-X . P  , 

with  y(0,0)  =  (x*,u*,w*)  ,  and  such  that  for  any  (e,r)  near  (0,0) 
and  r  >  0  ,  x(t,r)  is  a  locally  unique  unconstrained  local  minimizing  point 
of  W(x,t,r)  ,  with  gi[x(e,r),£]  >  0  ,  i-l,...,m  .and  v2w[x(e,r),£,r] 
positive  definite. 

The  relevance  of  equations  (2.1)  is  the  fact  that,  under  the  given 
conditions,  when  r  *  0  ,  they  are  necessary  conditions  that  must  hold  at  a 
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local  solution  of  P(0)  and,  with  r  >  0  ,  they  are  necessary  conditions 
for  an  unconstrained  minimum  of  W(x,f.,r)  .  The  latter  fact  can  be  made 
obvious  by  solving  for  u^  and  w^  in  (2.1)  and  obtaining 

V  L(x,u,w,e)  H  V  f  -  E  u  V  g  +  £  w.V  h 
X  x  i  x  i  j  X  j 

-  Vxf  “  r  I(l/gi)  Vxgi  +  (1/r)  E  h.V^ 

=  VxW(x,e,r)  . 

Thus,  if  y(e,r)  is  a  solution  of  (2.1),  then 

VxW[x(e,r),t,r]  «  VxL[x(e,r) ,u(e,r)  ,w(f.,r) ,e]  *  0  .  (2.2) 

This  explicit  connection  between  the  optimality  conditions  of  local 
solutions  of  P(e)  and  unconstrained  minima  of  W(x,e,r)  makes  it  possible 
to  approximate  information  characterizing  a  local  solution  of  P(e)  by 
algorithmic  calculations  associated  with  utilizing  W(x,e,r)  to  solve  P(e)  . 

In  particular,  differentiating  (2.2)  with  respect  to  e  yields 

V2^  x  +  V2  W  -  0 
x  e  ex 

and  using  the  fact  that  Vx*  is  positive  definite  (a  conclusion  of  Lemma  2.2) 
yields 

V  x  -  -V2W“V  W  (2.3) 

e  x  ex  v  ' 

evaluated  of  course  at  x(e,r)  .  Given  V^xU.r)  ,  the  derivatives  of  the 

multipliers,  Vu1(e,r)  and  Vcw^(e,r)  ,  can  then  be  calculated  by  differentiating 

the  last  two  systems  of  equations  of  (2.1)  at  x(e,r)  with  respect  to  e  . 

Lemma  2.3  (Fiacco  [10]).  For  e  in  a  neighborhood  of  e  *  0  , 
it  follows  that: 

(a)  lim  y(e,r)  -  y(e,0)  -  y(e)  ,  the  Kuhn-Tucker  triple 
r-K) 

characterized  in  conclusion  (b)  of  Lemma  2.1,  and 

(b)  lim.  V  y(e,r)  -  7  y(e,0)  -  V  y(e)  . 
r-K)  tee 
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Under  the  eondittons  of  Lemma  2.1,  x(i,r)  is  a  locally  unique  minimizing 
point  of  W(x,i  ,r)  .  Define  the  "optimal  value  penalty  function"  as 

W*(t  ,  r)  W(x(>  ,  r) , f , r)  . 

Armacost  and  Fiacco  [2,  Theorem  4  and  Corollary  4.1]  have  obtained  the  further 
results  useful  for  estimating  the  first-  and  second-order  sensitivity  of  the 
optimal  value  function  f*(f;)  of  Problem  P(c). 

Lemma  2 ■ 4  (Armacost  and  Fiacco  | 3 ] )  If  the  assumptions  A^  -  A^ 

fiold  for  Problem  P(<),  then  in  a  neighborhood  of  >  =  0  , 


(a) 

lira 
r  K1 

W*(f  ,r)  =  f*(t) 

(b) 

V  w* 
c 

(:  ,r)  =  Vt  L(y,^  ) 

(c) 

lim 

rK) 

V  W*(<:,r)  =  V  f 

(d) 

V^W* ( i" ,  r )  *  Vf2L(y(f, 

(e) 

lim 

r+0 

V2W*(f- ,  r)  =  V2f 

and 


wiiere  convergence  is  component  by  component  in  all  cases. 

Lemmas  2.1  -  2.4  enable  us  to  calculate  an  estimate  of  y(t),  Vfy(t)  , 

V  f*(t.)  ,  and  when  r.  is  near  0  and  r  is  near  0,  once  y(>:,r) 

is  available. 

In  the  next  section  we  briefly  present  the  algorithmic  implementation 
of  some  of  the  above  results. 

3.  THE  ALGORITHM 

The  penalty  function  algorithm  SUMT  estimates  the  solution  of  the  general 
mathematical  problem  P(e)  by  estimating  the  unconstrained  minima  of  the  penalty 
function  W(x,c,r)  at  successively  decreasing  values  of  the  penalty  function 
parameter  r  >  0  .  Under  conditions  weaker  than  those  assumed  here,  Fiacco 
and  McCormick  [13]  have  shown  that  as  r  approaches  zero,  the  sequence  of  the 
unconstrained  minima  of  W(x,c,r)  will  approach  a  solution  of  P(r)  .  Each 
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unconstrained  penalty  function  minimization  is  thus  a  "subproblem"  associated 
with  a  particular  value  of  the  penalty  function  parameter  r  . 

The  successive  steps  of  the  algorithm  for  sensitivity  analysis  that  are 
presented  in  reference  (1]  are  listed  below.  The  notation  x  or  x(t )  denotes 
the  estimate  of  a  solution  point  of  P(c)  calculated  by  SlfMT  for  a  given  value 
of  the  penalty  function  parameter  r  ,  where  £  denotes  the  value  of  the  problem 
parameter  for  which  this  sensitivity  is  estimated. 

Step  1:  Compute  a  representation  of  V^W  ^  =  V2W(x,e,r)  ^  by 

L-U  decomposition  using  the  SUMT  subroutines.  If  V2W  is  not 

positive  definite,  terminate  the  sensitivity  analysis. 

T 

Step  2:  Estimate  3(V  W  )/3e.  using  the  central  differencing 

*  3 

formula 

3(V  WT)/3£.  »  (1/2A)  (V  W(x,e+Ae.  ,r)T  -  V  W(x,c-Ve  ,r)T)  . 

x  j  x  j  x  j 


(Alternately,  the  matrix  of  partial  derivatives  could  be  determined 
analytically  from  W(x, E,r)  .) 


Step  3: 


Step  4: 


Step  5: 
using 


Calculate 

3 x(c)/3e.  *  -vV13('7  WT)/3e.  . 

j  x  x  j 

Estimate  V^g^^Cx.e)  and  V£,hi(x,e)  using 
3gi(x,e)/3e^  -  <l/2A)(gi(x,e+Ae^)  -  gi(x,e-Ae^))  ,  and 
3h1(x,£)/3ej  »  (l/2A)(hi(x,e+Aej)  -  h^x.E-Aej))  . 
Estimate  the  components  of  V£u(e)  for  i  »  l,...,m 


3ui(e)/3£j  -  -(r/gi(x,c)2)(Vjtg1(x,r)3x(£)/3£j 


Step  6;  Estimate  the  components  of  V£w(c)  for  i  ■  l,...,p 
using 

-  (l/r)(Vxh1(x,£) 3x(£)/3r.j 
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There  are  two  methods  for  estimating  V  f*(t.)  ;  the  first  using 
V^f *  =  f  V(  x  +  V(  f  ,  with  V^x  obtained  from  Step  3,  and  the  second 

method  using  the  gradient  of  the  penalty  function  W  ,  or  equivalently, 
the  Lagrangian  taken  with  respect  to  the  parameters  (Lemma  2.1,  conclusion 
(d).  Lemma  2.4,  conclusion  (b)).  Both  are  incorporated  in  the  computer  pro¬ 
gram  but  used  for  different  purposes.  The  former  method  gives  the  most 
accurate  estimate  of  V  f *(c )  and  is  summarized  in  Steps  7  and  8. 

Step  7:  Estimate  the  components  of  Vif(x,t)  using  the 
central  differencing  formula 

9f(x,t)/3e.  =  (1/2A)  (f (x.e+Ae^)  -  f(x,t>Ac.)  . 

Step  8:  Calculate  an  estimate  of  the  components  of  V  f*(c) 
using  the  results  of  Steps  3  and  7  as 

3f*(i)/3e.  =  V'  f(x,e)9x(t:)/3f_.  +  3f  (x,l) /3t  .  . 

J  x  3  J 

The  second  method,  using  the  gradient  of  the  Lagrangian  to  estimate 
V.f*(f),  is  computationally  less  expensive  and  is  used  to  obtain  rough  estimates 
that  single  out  the  more  crucial  parameters  for  further  analysis.  This  approxi¬ 
mation  is  calculated  as  follows. 

Step  9.  Estimate  the  components  of  V(_f*(c)  using  the  results  of 
Steps  4  and  7  as 

9f*(e)3Cj  =  3f(x,t)/3e.. 

m 

-  I  u,  (f.,r)3g,  (x,£)/3k 
i=l  1  1  J 

P 

+  Z  wi(t,r)3h  (x,e)/3c  . 

i=l  J 

4.  PROBLEM  DESCRIPTION 

The  nonlinear  programming  model  for  the  optimization  of  the  design  of 
a  vertically  corrugated  transverse  bulkhead  of  an  oil  tanker,  discussed  in 
Chapter  6  of  [8],  was  selected  for  the  sensitivity  analysis.  The  material  in 
this  section  follows  the  presentation  in  [8]  rather  closely. 
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Vertical  transverse  bulkheads  form  the  lateral  walls  of  the  internal 
compartments  of  the  tankers  that  hold  liquid  cargo.  Figure  1  shows  three 
different  views  of  a  corrugated  bulkhead,  consisting  of  a  bottom,  middle, 
and  top  panel,  fastened  together  along  the  stringers  EF  and  CD.  It  is  assumed 
that  the  shape  of  the  corrugations  are  identical  in  all  panels  but  their  thick¬ 
nesses  are  allowed  to  vary  from  panel  to  panel.  The  lengths  of  the  top,  middle, 
and  bottom  panels  denoted  by  £  ,  £  ,  and  ,  respectively,  are  fixed,  as 

is  the  common  width  B  of  the  panels. 

The  problem  design  variables  are  indicated  in  Figure  2  and  defined  as 
follows : 

b^  -  width  of  the  flange 
*  length  of  the  web 

d  *  depth  of  corrugation 

tt  ■  thickness  of  top  panel 

t  *  thickness  of  middle  panel 

TR 

tfa  -  thickness  of  bottom  panel. 


The  width  of  a  corrugation  s  is  depicted  and  used  in  the  model  for 

convenience.  It  is  not  a  design  variable,  being  determined  once  the  indicated 

2  2  1/2 

variables  are  specified.  In  fact,  it  is  easy  to  see  that  s  =  b^  +  ^ 

The  design  parameters  involved  in  the  formulation  of  the  model  and  their 
specified  values  are  listed  in  Table  1.  Those  characterizing  the  shape  of  the 
bulkhead  are  depicted  in  Figures  1  and  3.  Figure  3  also  indicates  the  load 
configuration  which  will  enter  into  the  constraints  of  the  model  to  be  developed 
subsequently. 

Objective  Function 

It  is  desired  to  determine  values  for  these  design  variables  that 
minimize  the  weight  of  the  corrugated  bulkhead  subject  to  satisfying  a  number 
of  constraints  to  be  considered  subsequently.  The  total  weight  w  of  the 
panels  is  easily  calculated  to  be 


Figure  1.  Vertical  corrugated  transverse  bulkhead. 
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Table  1 

PROBLEM  DES1CN  PARAMETERS 
(MODEL  INPUT  DATA) 


HI 

Parameter 

Meaning 

Value 

n 

r 

Weight  per  unit  volume  of  the  material 

7.85ton/cm3 

B 

Width  of  the  panel 

476  cm 

B 

l 

t 

Length  of  the  top  panel 

495  cm 

H 

i 

m 

Length  of  the  middle  panel 

385  cm 

5 

h 

Length  of  the  bottom  panel 

315  cm 

6 

h 

a 

Distance  between  free  liquid  level  and 
top  of  the  structure 

250  cm 

7 

h 

t 

Distance  between  free  liquid  level  and 
middle  of  the  top  panel 

498  cm 

8 

h 

m 

Distance  between  free  liquid  level  and 
middle  of  the  middle  panel 

938  cm 

9 

\ 

Distance  between  free  liquid  level  and  middle 
of  the  bottom  panel 

1288  cm 

10 

hu 

Distance  between  free  liquid  level  and  base 
of  the  top  panel 

745  cm 

11 

“lm 

Distance  between  free  liquid  level  and  base 
of  the  middle  panel 

1130  cm 

12 

hlb 

Distance  between  free  liquid  level  and  base 
of  the  bottom  panel 

1443  cm 

i 

13 

min 

ct 

Minimum  allowable  thickness  of  the  top  panel 

1.05  cm  j 

14 

tmin 

m 

Minimum  allowable  thickness  of  the  middle  panel 

1 

1.05  cn 

15 

min 

b 

Minimum  allowable  thickness  of  the  bottom  panel 

1.05  cm 

16 

e 

Effectiveness  of  the  flange  (dimensionless) 

0.8 

17 

ki 

Coefficient  (function  of  maximum  allowable 
bending  stress) 

6.94  x  10  ^cm  ^ 

18 

k2 

Corrosion  coefficient 

0.15  cm 
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t  i+t  H  +t,  l, 

w  =  rB(b.+b_)  — ^--m— (1) 

12s 

where  P  is  the  weight  per  unit  volume  of  the  material. 

It  is  assumed  that  the  design  variables  are  to  be  given  in  centimeters. 
I‘  will  be  given  in  tons  per  cubic  centimeter. 


Geometrical  Constraint 

As  indicated  in  Figure  2,  geometrical  limitations  suggest  that  the  length 
of  the  web  should  at  least  be  equal  to  the  depth  of  the  corrugation.  Therefore, 
it  is  required  that 

b2  -  d  _>  0  (2) 


Bending  Stress  Constraint 


Consider  any  panel  of 
distance  from  load  level  zero 
stress  at  the  supports  of  any 
given  by 


a 


the  type  described  and  denote  by  h  ,  the  vertical 
to  the  middle  of  the  given  panel.  The  bending 
such  panel  of  thickness  t  and  length  l  is 


(kg/cm  ) 


where  e  (the  effectiveness  of  the  flange)  is  a  dimensionless  constant  and 

3 

Y  =  .001  kg/cm  is  the  specific  gravity  of  fresh  water.  (The  numerator  in  the 
expression  for  a  is  the  bending  moment  at  the  supports  in  kg  cm  ,  and  the 
denominator  is  the  section  modulus  in  cm^  .) 

The  maximum  permitted  bending  stress  is  given  to  be  1200  kg/cm2  , 
thus  requiring  that 

a  <  1200  (kg/cm2)  . 


Therefore,  the  prior  expression  for  o  and  the  above  inequality  result  in 
the  constraint 
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where  k^ 


td 

2 


+  ble) 


k^hsJi 


.694  x  10  7  (cm 


(3) 


Moment  of  Inertia  Constraint 


An  additional  constraint  on  the  given  panel  is  imposed  as  follows: 


(!)  (¥)  ft  *  v)  > 


2.2 


(^hsJL2)4^ 


(4) 


where  the  expression  on  the  left  is  the  moment  of  inertia  of  the  panel. 


Thickness  Requirement  Constraint 

For  a  given  panel  it  is  required  that 

1tmin 

(.39) (1.05  bj)  /.01  hL  +  k?  (5) 

(.39)  (1.05  b2)  /ToTh^V  k2 

where 

t  =  plate  thickness  (centimeters) 

tmin  *  minimum  allowable  plate  thickness  (centimeters) , 
a  function  of  ship  length 

h^  ■  height  to  load  level  zero  from  the  bottom  of  the 
panel  centimeters,  and 

k2  “  corrosion  allowance  (centimeters). 

Hatural  Constraints 

Since  the  design  variables  b,,  b„,  d,  t ,  t  and  t.  are  measures  of 

i  t  t  m  t> 

length  (dimensions  of  the  corrugated  bulkhead),  they  are  constrained  to  be 
nonnegative. 
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Design  parameters  h,  £,  tmln  and  in  constraints  (3),  (A),  and 

(5)  are  subscripted  by  t,  m  and  b  to  define  the  relative  constraints  for 
the  top,  middle,  and  bottom  panels,  respectively. 

The  Design  Problem 

Introducing  variables  x^  through  x^  for  the  design  variables 
bj,  b^,  d,  t^,  t^,  tj  ,  respectively,  and  denoting  the  objective  function  w 
by  f  ,  yields  the  following  nonlinear  programming  problem. 


minimize  f  rB(xn+x  )  (£  x.+R.  x  +£  x. )  [x  +  (x2-x2)2]  1 
Jz  t  4  m  o  d  o  L  2  3 


subject  to 


x2x3x4  +  3e  X1X3XA  ~  6klht£t^Xl  +  (x2-x2) 2 ]  >  0 


2  2. 


13  5  1  m  m  1  T  VA2  A3' 


2  3  6 

2 


’7  ‘  2  3  6 


C1X3X6 


min 


Constraint 

Number 


x2  "  x3  >  0  •  (1) 


X2X3X5  +  3e  xi x,tx  s  ~  6k1hm^fx1  +  x-a ) 2 J  >  0  • 


(2) 

(3) 


3e  x^x^  -  +  (x^-x2^] 

>  0  . 

(4) 

26.4(k1ht^)4/3[x1  +  (x2-x2)I]4/3 

>  0  . 

(5) 

26'4(klVm>4/3fxl  +  (x2-x3^14/3 

>  0  . 

(6) 

26*4(klIb£b)4/3[xl  +  (x2"x3>^4/3 

>  0  . 

(7) 

$ 
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x4  -  [.39  •  1.05(.01hlt)2(.01x1)  +  k2l  >  0  .  (9) 


no 


x4  -  [.39  •  1.05(.01hlt)2(.01x2)  +  k2]  >  0  .  (10) 


’ll 


min 

*5  ~  tm  ^  0  •  <U> 


’12 


’13 


x5  -  [.39  •  1.05(.01hlm)2(.01x1)  +  k9]  >  0  .  (12) 

1 

x5  -  [.39  •  1.05(.01hlm)2(.01x  )  +  k  ]  >  0  .  (13) 


*14 


x6  ~  t™in  >  0  .  (14) 


*15 


’16 


X6  -  [.39  ♦  1.05(.01hlb)2(.01x1)  +  k2]  >  0  .  (15) 


x6  -  [.39  •  1.05(.01hlb)  (.01x2)  +  k2J  >  0  .  (16) 


817  "  g18 


x.  >  0  ;  1  =  1  and  3  .  (17-18) 


5.  COMPUTATIONAL  RESULTS 

Phase  1:  Solution  of  the  Problem  P(e) 

This  phase  involves  the  solution  of  the  problem  without  parameter 
perturbation.  Denote  by  e  the  vector  of  desigu  parameters  listed  in 
Table  1,  i.e.,  the  ith  component  of  c  designating  the  ith  parameter  in 
Table  1.  With  this  notation,  let  c  represent  the  vector  whose  components 
are  the  respective  data  given  in  Table  1.  Therefore,  the  design  problem 
given  in  the  previous  section,  with  the  problem  design  parameters  equal 
to  the  values  given  in  Table  1,  is  denoted  P(c),  in  conformance  with  the 
notation  of  Section  2. 

Table  2  gives  the  calculated  solution  x(e)  (to  two  significant 
figures)  of  P(c)  ,  along  with  the  solution  presented  in  [8].  Our  optimal 
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TABLK  2 


OPTIMAL  SOLUTIONS  FOR 
PROBLEM  P(0 


Starting 

Point 

(cm) 

Solution  in  [ 7 j 
(cm) 

Our 

Solution 
x(t)  (cm) 

xr  =  bL  =  45.30 

57.80 

57.82 

x2  =  b9  =  43.20 

57.80 

57.82 

x^  =  d  =  30.50 

37.80 

35.69 

x4=tt=  1.20 

1.05 

1.05 

xc  =  t  =  1.20 

5  m 

1.05 

1.05 

x  =  r  =  1.30 

o  b 

1.05 

1.05 

f  =  w  =  6.40  tons 

5.34  tons 

5.25  tons 
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solution  matches  that  of  reference  [8]  except  for  x^  (the  corrugation  depth). 

This  slight  difference  in  the  value  of  ,  and  consequently  in  the  value  of 

objective  function  W  is  believed  to  be  due  solely  to  the  difference  in 
manipulation  of  the  input  data  to  the  model.  In  the  current  study  the  design 
parameter  values  were  manipulated  separately  in  the  model,  as  encountered  in 
a  given  calculation.  The  results  given  in  [8]  probably  involved  the  calcula¬ 
tion  of  coefficients  (which  are  themselves  functions  of  design  parameters) 
"externally,"  using  the  aggregate  results  as  input  data  to  the  model.  The 
rounding  of  coefficients  so  calculated  could  readily  explain  the  resulting 
minor  discrepancy. 

Phase  2:  Sensitivity  Analysis  of  Problem  P(£)  for 
Right  Hand  Side  Constraint  Perturbations 

The  problem  P(e)  with  right  hand  side  perturbation  a  of  the  constraints 
will  be  called  Problem  R(a) .  Obviously,  R(0)  =  P(e),  so  the  optimal  solution 
x(t;)  of  P(e)  given  in  Table  2  is  also  x(0)  ,  the  optimal  solution  of  R(0)  . 

Denote  by  a  =  c  the  vector  a  ,  all  of  whose  components  are  equal 
to  the  number  c  .  It  was  desired  to  analyze  the  solution  of  R(a)  for 
cx  =  +  .25  and  a  =  +  .50  . 

The  optimal  solutions  x(a)  and  the  first  order  sensitivities 
V(jf*(a)  ,  V(<x(a)  and  V^u(a)  were  obtained  for  the  given  values  of  a  . 

The  following  results  were  also  obtained. 

(1)  The  assumptions  of  Lemma  2.1  were  shown  to  hold  at  x(0) 

and  the  results  of  the  Lemma  were  computationally  corroborated, 
for  a  ■  +  .25  and  a  ■  .50  . 

(2)  Using  the  solution  and  first  order  sensitivity  information 
associated  with  Problem  R(0) ,  the  optimal  solution  for  the 
perturbed  problems  R(+.25)  and  R(+.50)  were  estimated 
and  compared  to  those  obtained  by  direct  solution. 

The  following  elaborates  on  the  above  results  (1)  and  (2).  Table  3 
depicts  the  Lagrange  multipliers,  constraint  values  and  sensitivity  informa¬ 
tion  associated  with  the  optimal  value  function.  It  is  clear  from  this  table 
that  x(0)  is  feasible,  the  multipliers  are  nonnegative  and  complimentary 
slackness  is  satisfied  (to  two  significant  places). 
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TABLE  3 


LAGRANGE  MULTIPLIERS  AND  CONSTRAINT 
VALUES  FOR  PROBLEM  R(0)  AT  x(0) 

(Excluding  non-negativity  constraints) 


(i) 

Constraint 

No. 

Value  of 
Constraint  g^ 

Value  of 
Lagrange 
Multiplier  u^ 

3f*(0) 

Da. 

i 

-2 

-2 

1 

22.13 

.4518784  x  10 

.4518777  x  10 

-4 

-4 

2 

2117.22 

.4723479  x  10 

.4723564  x  10 

-4 

-4 

3 

1385.41 

.7218857  x  10 

.7218987  x  10 

-4 

-4 

4 

1868.92 

.5351946  x  10 

.535204  x  10 

5 

41993.15 

.2381633  x  10“5 

.2381608  x  lO-^ 

*  6 

.043 

1 

.2300013  x  10 

.2300013  x  10l 

-5 

-5 

7 

27945.57 

.3580375  x  10 

.3580337  x  10 

*  8 

.48  x  10-5 

.2078305  x  107 

.2078305  x  107 

9 

.25 

.3940568 

.3940568 

10 

.25 

.3940570 

.3940570 

*  11 

.96  x  10-5 

.1040634  x  107 

.1040634  x  107 

12 

.10 

.9604345 

.9604345 

13 

.10 

.9604360 

.9604360 

*  14 

.49  x  10'4 

.2029215  x  106 

.2029215  x  106 

*  15 

.31  x  10"4 

.3221582  x  106 

.3221582  x  106 

-4 

6 

6 

*  16 

.13  x  10 

.7974950  x  10 

.7974950  x  10 

Constraints  marked  by  asterisk  (*)  are  binding. 
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Based  on  Table  3  it  was  concluded  that  the  constraints  g, ,  g„,  g,  , 

o  o  11 

814*  815  anc*  816  were  binding.  Rather  large  value  of  the  constraint  g^ 

compared  to  the  values  of  the  remaining  binding  constraints  may  encourage  the 
skeptical  reader  to  disagree  with  the  fact  that  this  constraint  is  binding. 
However,  one  should  notice  that  the  coefficients  involved  in  constraint  g^ 

are  very  large  compared  to  the  coefficients  of  the  remaining  binding  con¬ 
straints  (this  is  also  reflected  in  the  values  of  the  components  of  Vg^ 

listed  in  Table  4).  By  proper  scaling  of  this  constraint  one  may  reduce  its 
value  to  the  order  of  the  values  of  the  remaining  binding  constraints.  A 
sufficient  reason  that  constraint  g^  may  be  concluded  to  be  binding  is  that 
the  optimal  value  function  is  sensitive  to  parameter  (to  be  shown  in 

Phase  3)  while  k^  appears  neither  in  the  objective  function  nor  in  any  of 

the  binding  constraints  but  constraint  g,  .  Incidentally,  this  gives  an 

o 

interesting  example  of  the  sort  of  insight  that  is  provided  by  a  thorough 
sensitivity  analysis. 

Table  4  depicts  the  gradients  of  these  binding  constraints  and  their 
corresponding  optimal  Lagrange  multipliers  at  the  solution  point  for  Problem 
R(0).  It  is  clear  that  the  binding  constraint  gradients  are  linearly  inde¬ 
pendent  (A3)  and  the  associated  Lagrange  multipliers  are  positive.  Therefore 
strict  complimentarity  slackness  (A^)  holds  at  x(0)  .  Since  x(0)  was 
calculated  to  be  a  stationary  point  of  the  Lagrangian  function,  it  follows 
from  these  facts  that  (A2)  is  satisfied.  This  also  means  that  x(0)  is 
indeed  an  isolated  local  solution.  Finally  the  problem  functions  are  all 
twice  differentiable  (A^).  Therefore,  relative  to  two  significant  figures 
of  accuracy,  the  four  assumptions  of  the  Lemma  2.1  have  been  verified  for 
Problem  R(0).  So  the  results  of  that  lemma  hold  for  Problem  R(a) ,  with  a 
near  zero. 

The  last  column  in  Table  3  depicts  the  sensitivity  of  the  optimal 
value  function  for  Problem  R(0) .  As  it  is  seen  the  entries  in  this  column 
closely  correspond  to  the  respective  Lagrange  multipliers,  listed  in  column  3 
of  Table  3,  as  expected  from  result  (d)  of  Lemma  2.1. 
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TABLE  4 


GRADIENTS  OF  THE  BINDING  CONSTRAINTS 
AND  CORRESPONDING  LAGRANGE  MULTIPLIERS 
FOR  PROBLEM  R(a)  AT  x(a) 


V*6 

V% 

V*ll 

V«14 

V*15 

'816 

-288 

0 

0 

-.01557 

0 

-3352 

0 

0 

-.01557 

r(-.T>) 

15752 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

313158 

0 

1 

0 

0 

0 

0 

0 

0 

l 

1 

1 

u(~. 25) 

2.53 

2.16  x  106 

1.98  x  105 

3.23  x  105 

8.51  x  105 

-183 

0 

0 

0 

-.01557 

0 

-2975 

0 

0 

0 

0 

-.01557 

K(O) 

17393 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

250354 

0 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

u(O) 

2.30 

2.07  x  106 

1.04  x  106 

2.03  x  105 

3.22  x  105 
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Tabic  5  depicts  the  optimal  solutions,  optimal  values,  and  a  typical 
Lagrange  multiplier  ug(°0  for  problems  R(a)  ,  where  a  =  +.25,  +.50  , 

obtained  by  direct  calculation.  Figures  4,  5  and  6  suggest  that  x(«)  , 

f*(a)  and  uQ(a)  are  continuous  and  smooth  functions  of  c  for  the  entire 
o 

range  of  parameters  considered  (a  result  known  to  hold  for  n  near  0  ,  by 

T 

conclusion  (b)  of  Lemma  2.1),  where  a  =  c  =  c(l,1,...,l) 

Table  6  indicates  the  binding  constraints  (using  two  significant 
figures)  for  Problem  R(a)  at  x(a)  .  As  may  be  seen,  the  binding  constraints 
for  all  perturbed  problems,  except  R(-.50)  ,  remain  the  same.  Moreover, 

Table  4  shows  that  the  gradients  of  these  binding  constraints  are  linearly 
independent  and  we  determined  that  strict  complimentarity  slackness  holds  true. 
So,  the  results  of  Lemma  2.1  must  also  hold  for  the  problems  R(+.23)  and 
R(+.50)  .  The  change  in  the  binding  constraint  index  set  for  problem  R(-.50) 
means  that  the  perturbation  a  =  -.50  is  not  small  enough  to  retain  the  general 
solution  structure  of  the  problem. 

In  order  to  estimate  the  optimal  solution  of  Problem  K(a)  based  on 
information  from  Problem  R(0)  the  following  first  order  (Taylor's  Series) 
extrapolation  formula  was  used, 

y(a)  =  y(0)  +  V^y(0)  •  a 

where  y(0)  and  Vay(0)  are  the  estimates  of  a  Kuhn-Tucker  Lriple  and  its 
derivatives  with  respect  to  a  for  problem  R(0) .  Although  the  magnitude  of 
the  perturbations  were  quite  significant,  the  extrapolated  values  are  in  close 
agreement  with  those  of  direct  calculation.  Extrapolation  was  also  done  for 
Problem  R(-.50)  to  study  the  discrepancy  between  extrapolated  and  calculated 
results.  Table  7  shows  these  results  and  Figure  7  depicts  some  of  these, 
namely  x^(a)  and  x^(ot)  ,  graphically. 

In  order  to  computationally  check  the  stability  of  the  algorithm's 
estimates  of  solution  and  sensitivity  information  of  Problem  R(ct),  f(*),  x(*) 
u(»)  and  their  partial  derivatives  with  respect  to  a  were  recorded  along 
the  algorithm  solution  trajectory.  Values  along  the  trajectory  path  are  plotted 
In  Figure  8  for  a  sample  of  typical  values,  for  the  Problem  R(.50).  As  it  is 
seen,  the  solution  estimates  are  virtually  constant  from  the  third  subproblem 
onwards . 
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•p\  -to 


Figure  4.  x(a)  versus  c  for 

•  for  x^  and 
0  for  x^ 

x  for  x.and  x.  Common 
4  5 


Problem  R(ot)>  where  a  =  c  . 

+  for  x, 
o 

*  for  x^,  Xj  and  x^  Common 
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TABLE  6 


BINDING  CONSTRAINTS  FOR  PROBLEM 
R(oO  AT  x(a) 


Constraint^. 

-.50 

-.25 

0 

+  .25 

+  .5( 

81 

- 

- 

- 

S2 

- 

- 

g3 

B 

- 

- 

- 

- 

g4 

B 

- 

- 

85 

- 

- 

g6 

B 

B 

B 

B 

g7 

- 

- 

- 

- 

g8 

B 

B 

B 

B 

B 

g9 

- 

- 

- 

- 

gL0 

- 

- 

- 

8ll 

- 

B 

B 

B 

B 

g12 

- 

- 

- 

- 

- 

g13 

- 

- 

- 

814 

- 

B 

B 

B 

B 

g15 

B 

B 

B 

B 

B 

816 

B 

B 

B 

- - -  ■ 

B 

B 

B  stands  for  binding  constraints 
-  stands  for  nonbinding  constraints 


where 
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TABLE  7 


EXTRAPOLATED  SOLUTIONS  VERSUS  CALCULATED  SOLUTIONS 
FOR  PROBLEM  R(a) 


R(~. 50) 

R(-.25) 

R(0) 

R(. 25) 

R(. 50) 

e1 

x!  ^ 

57.82 

57.82 

- 

57.82 

57.82 

58.45 

57.82 

57.82 

57.82 

57.82 

e 

57.82 

57.82 

- 

57.82 

57.82 

58.45 

57.82 

57.82 

57.82 

57.82 

e 

*3  c 

42.88 

39.28 

- 

28.49 

39.91 

35.69 

1 

30.10 

e 

.55 

.30 

- 

1.30 

1.55 

.55 

.80 

1.05 

1.30 

1.55 

H 

.55 

.80 

- 

1.30 

1.55 

.61 

.80 

1.05 

1.30 

1.55 

.55 

.80 

- 

1.30 

1.55 

.56 

.80 

1.05 

1.30 

1.55 

e 

f* 

c 

1 

4.137 

- 

6.358 

7.469 

B 

4.145 

5.248 

6.356 

7.467 

^e,  extrapolated 
2 

c,  calculated 
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6  Calculated 
+  i^xtrupolat^d 
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Numerous  inferences  can  be  drawn  from  such  information.  For  example, 
note  that  the  constraints  gg,  g  ^  and  g^  whicti  correspond  to  the  minimum 

plate  thickness  requirement  for  top,  middle,  and  bottom  panels  remain  binding 
for  the  perturbed  problems  over  the  range  of  perturbations  considered  (i.e., 

-.25  to  +.50).  On  the  other  hand,  because  the  weight  of  the  bulkhead  is 
linearly  proportional  to  the  plate  thickness,  any  relaxation  on  minimum  thick¬ 
ness  requirement  would  decrease  the  bulkhead  weight  (objective  function)  by  a 
considerable  amount.  Based  on  the  results  given  in  Table  3,  the  saving  on 
bulkhead  material  weight  will  be  about 

.207  +  .104  +  .020  =  .331  Lons 

per  one  millimeter  relaxation  of  the  required  plate  thickness.  Figure  5,  which 
shows  the  optimal  function  values  versus  perturbations,  confirms  this  observation. 


Phase  3:  Sensitivity  Analysis  with  Respect  to 
Design  Parameter  Perturbations 

This  study  phase  seeks  to  obtain  sensitivity  information  when  the  design 
parameters  are  perturbed.  As  noted.  Table  1  lists  all  the  design  parameters 
and  their  initial  values. 

The  parameters  1  (density  of  steel),  e  (dimensionless  constant  in 
section  modulus  formula),  B  (the  width  of  the  panel)  and  ha  (the  distance 
from  top  of  the  bulkhead  to  free  liquid  level)  are  assumed  to  be  subject  to 
negligible  variation  and  so  their  perturbations  are  excluded  from  the 
sensitivity  analysis. 

It  is  important  to  note  that  the  design  parameters  £.  ,  %  ,  f,  ,  h  ,  h 

t  m  b  i  m 

*'b’  ''it*  ^lm  an<*  ^lb  3re  not  in^ePendent.  for  example,  as  can  easily  be 

verified  from  Figure  3,  these  parameters  may  be  expressed  in  terms  of 

I  .,  i  and  h,,  as  follows: 
t  w  1b 


h1K  -  (250  +  £  +  i  ) 
lb  t  m 


hlt  =  250  +  f 


h  =  250  +  l  +  ^2- 
m  t  2 


hb  "  125  + 


H„+IL  +h.. 
t  m  lb 
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a. 


i. 


k 


hlt  =  250  +  Jit 

h.  =  250  +  £  +  Jl 
lm  t  m 


In  this  analysis  of  the  problem,  we  explicitly  impose  these  equations 
and  express  all  the  design  parameters  as  functions  of  and  h^  »  as 

indicated.  This  eliminates  unnecessary  sensitivity  calculations  and  strictly 
accommodates  the  dependencies. 


The  "independent  parameters" 
shown  in  Table  8.  For  convenience. 


on  which  the  sensitivity  is  analyzed  are 
we  relabel  the  jth  parameter  c.  . 


As  previously  noted,  the  corresponding  numerical  values  of  e.  are  given  in 

Table  1.  These  values  are  taken  as  the  initial  value  of  the  problem  parameters 
and  are  again  designated  by  r.  . 


Aside  from  an  analysis  of  the  Problem  P(e),  the  problem  was  also  analyzed 
under  six  different  perturbations.  Namely,  all  of  the  independent  design 
parameters  were  perturbed  simultaneously  by  the  same  fraction  of  their  original 
values.  The  perturbations  were  +1,  +2  and  +5  percent  of  the  original  problem 
parameter  vector  t:  .  The  corresponding  problems  are  hence  designated  by 
P(l.01e),  P(,99e),  P(1.102e),  P(.98e),  P(1.05e)  and  P(.95e),  respectively, 
following  our  usual  notation. 


After  calculating  the  optimal  solutions  and  corresponding  first  order 
sensitivity  information  for  the  perturbed  problems,  the  following  analysis 
was  completed: 

(1)  For  small  perturbations  the  results  of  Lemma  2.1 
for  Problem  P(e)  were  verified, 

(2)  Using  the  solution  and  the  first  order  sensitivity 
information  of  the  Problem  P(e),  optimal  solutions 
for  the  Problem  P(c)  were  estimated  and  compared 
with  those  obtained  by  direct  solution. 


Tables  8,  9,  and  10  show  the  sensitivities  of  the  optimal  value  function, 
solution  point  and  Lagrange  multipliers  with  respect  to  the  independent  design 
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TABLE  8 


SENSITIVITY  OK  f*(e)  WITH  RESPECT  TO  THE  INDEPENDENT 
DESIGN  PARAMETERS  e  FOR  PROBLEM  P (~ ) 


Parameter 

Index 

j 

Parameter 

c . 

J 

Sf*(e) 

i 

i 

t 

4393 

2 

m 

8582 

3 

htb 

349 

4 

min 

t 

2078305 

■j 

^min 

m 

1040634 

6 

min 

b 

20292 L 

7 

k2 

1119656 

TABLE  9 


SOLUTION  POINT  SENSITIVITY  Dx  /3e 
Fur  PROBLEM  P(t) 


.2  x  10  J 
.015842 
.020053 
.774  x  10” 
.1297  x  10 
64.2357 
64.23520 
.4453  x  10 


2 

3 

.2  x  10"5 

.  * 

-.013829 

•  10226 

-.020053 

-.3625  x  10"2 

726  x  10“3 

.119  x  10-3 

-.123  x  10"2 

-14.39778 

64.23601 

11.61322 

-64.23554 

-11.61314 

-.3887  x  10+8 

.2874  x  10+9 

5 

6 

- 

- 

- 

-.208  x  10~3 

- 

“A 

.11  x  10 

-4 

.999999 

-.18  x  10 

- 

.999925 

-4 

— 

.82  x  10  H 

-1831.179 

-.5848  x  10+6 

1 

*Sensitivities  with  less  than  10  in  absolute  value  are  indicated  by 
(-)  in  the  above  tablq* 


TABLE  10 

LAGRANGE  MULTIPLIERS  SENSITIVITY 
FOR  PROBLEM  P(e) 
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I'. unmet  rin  <>l  I'lulilrm  l’(  •  > .  An  shown  in  Table  M  l  lie  up  (.Inin  I  value  lieutlon 

is  extremely  sensitive  to  minimum  allowable  panel  thicknesses  tmin,  tmin,  tmin 

t  m  b 

and  corrosion  allowance  and  apparently  sensitive  to  the  allowable  bending 

stress  factor  k^  .  It  is  only  partially  sensitive  to  panel  lengths  #,  and 

and  l  ,  and  practically  insensitive  to  the  total  depth  of  panel  h,,  . 

m  lb 

It  is  well-known  that  in  order  to  obtain  a  meaningful  interpretation 
of  the  implications  of  a  sensitivity  analysis,  the  information  must  be  care¬ 
fully  analyzed  in  the  perspective  of  the  application.  For  example,  the 
parameter  derivatives  of  the  objective  function  f*(0)  are  the  instantaneous 
rates  of  change.  A  rate  of  change  can  itself  change  significantly  away  from 
the  base  point,  thus  rendering  extrapolations  invalid  for  finite  parameter 
changes.  Thus,  further  analyses  such  as  those  reported  here  are  essential. 
Another  caution  must  be  directed  to  equating  the  rate  of  change  with  a  change 
"due  to  a  one  unit  change"  in  the  given  parameter  value.  It  is  a  change  per 
unit  change,  but  only  at  a  given  parameter  value. 

A  unit  change  in  a  parameter  value  may  not  be  feasible  or  meaningful. 

For  example,  the  initial  value  of  the  parameter  is  .694  x  10~^  cm~l  .  A 

one  unit  decrease  is  not  possible.  A  "reasonable"  change  for  this  parameter 

-9 

might  be  expected  to  be  on  the  order  of  +10  .  Given  this  scaling,  it  follows 

^hat  the  objective  function  is  only  relatively  mildly  sensitive  to  Lhc  changes 
in  this  parameter.  Similarly  ,  the  sensitivity  of  th"  parameters  k2  and 
the  minimum  allowable  panel  thicknesses  should  all  probably  be  scaled 
by  a  factor  of  10  .  Appropriate  scaling  of  parameter  changes  depends  on 

the  context  of  the  application,  parameter  interactions  that  may  not  even  be 
explicitly  represented  in  the  model  and  other  considerations.  This  is 
ultimately  a  "management"  decision  and  is  beyond  the  scope  of  this  paper. 

We  shall  assume  that  specified  feasible  changes  are  stipulated  and  concen¬ 
trate  on  measuring  their  effects. 

Table  II  and  Figure  9  show  the  calculated  optimal  solutions  for  Problem 
P(r),  when  it  is  stipulated  that  the  given  parameter  values  be  altered  by  the 
indicated  fractions  of  their  original  values. 
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TABLE  11 

OPTIMAL  SOLUTIONS  FOR 
PROBLEM  P  (e) 


P(0.95e) 

P(0.98e) 

P(0.99e; 

— 

P(e) 

P(1.0lE) 

P(1.02c) 

P(l. 05c) 

(cm) 

63.02 

61.53 

57.25 

57.82 

58.42 

59.01 

61.64 

x2  (cm) 

63.04 

61.53 

57.26 

57.82 

58.43 

59.02 

61.65 

x3  (cm) 

38.93 

35.90 

35.19 

35.69 

36.2 

36.71 

38.41 

(cm) 

1.0 

1.03 

1.04 

1.05 

1.06 

1.07 

1.10 

x^  (cm) 

1.0 

1.03 

1.04 

1.05 

1.06 

1.07 

1.10 

x&  (cm) 

1.13 

1.10 

1.04 

1.05 

1.06 

1.07 

1.12 

u8<0 

1964924 

2009085 

2055230 

2078305 

2101457 

2125542 

2189718 

f*(t.)  (Tons) 

5.176 

5.175 

5.190 

5.248 

5.306 

5.365 

5.541 

.(A-.  KMHil 
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It  is  obvious  that  the  assumptions  of  the  main  sensitivity  theorem 
holds  true  for  P(£)  (because  P(f  )  is  identical  to  R(0)  and  we  have 
already  verified  the  validity  of  these  assumptions  for  R(0)  in  the  previous 
part).  So  the  results  of  this  theorem  should  also  hold  true  for  Problem  P(c) 
provided  that  the  introduced  perturbations  are  not  too  large. 

The  result  a  of  Lemma  2.1  is  immediate.  Figures  9  and  10  are  included 

to  illustrate  the  result  b  of  this  lemma.  As  depicted,  x(>)  and  u  (e) 

8 

(a  sample  of  optimal  Lagrange  multiplier)  appear  to  he  smooth  functions  of 

d  ,  for  fc!  between  0.99  and  1.01  ,  where  t  =  dr.  .  However  there  is  a 

definite  sharp  change  (or  discontinuity)  for  uc(t)  as  well  as  x  (i),  x,(f) 

»  i  l 

and  xfi(c)  for  s<jme  value  of  t  between  0.99c  and  0.  98>:  .  This  may 

be  due  to  the  fact  that  perturbations  more  than  -0.01c  (in  magnitude)  are 
not  sufficiently  small  to  keep  the  structure  of  the  perturbed  problem  es¬ 
sentially  unchanged.  Kxamining  the  binding  constraints  for  the  perturbed  prob¬ 
lems  confirms  this.  Table  12  shows  the  binding  constraints  of  Problem  P(( ) 
at  different  perturbations. 

As  shown  for  perturbations  below  -0.01c  .  The  binding  constraints  g, 

b 

and  g  are  well  satisfied  and  constraints  g^,  and  8^3  become  binding. 

Moreover  constraint  g^  is  no  longer  binding  for  perturbations  of  0.05c  . 

Furthermore,  as  depicted  in  Table  13,  the  gradients  of  binding  constraints  of 
Problem  P(c)  are  linearly  independent  and  strict  complimentarity  slackness 
holds  true  for  c  between  0.99c  and  1.02e  .  (Result  c  of  Lemma  2.1).  So 
all  the  results  of  Lemma  2.1  for  above  range  of  perturbation  (i.e.,  0.99c  to 
1.02c)  are  computationally  verified. 

As  before  estimates  were  made  for  the  optimal  solution  of  the?  perturbed 
Problem  P(e),  based  on  the  solution  and  sensitivity  information  of  the  unper¬ 
turbed  Problem  P(c),  by  first  order  extrapolation.  Table  14  and  Figure  11  show 
the  results  of  this  extrapolation  and  compare  them  with  the  results  of  direct 
calculation.  It  is  obvious  that  one  should  not  attempt  to  extrapolate  from  the 
solution  of  P(c)  to  solutions  of  P(  0.98e)  ,  P(0.95e)  and  P(1.05e)  .  We 
have  deliberately  calculated  this  extrapolation  to  show  the  consequences. 
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TABLE  12  '* 

BINDING  CONSTRAINTS  FOR  PROBLEM 
P(c)  AT  x(t) 
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2.8  x  10°  1.04  x  10"  2.03  x  10"  3.22  x  10"  7.97  x  10' 
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TABLE  13  (Cont'd) 

GRADIENTS  OF  THE  BINDING  CONSTRAINTS 
AND  CORRESPONDING  LAGRANGE  MULTI  I’Ll  ERS 
FOR  PROBLEM  P(e)  AT  x(e) 


V^6 

% 

**11 

V*!4 

%5 

V*16 

-194 

0 

0 

0 

-.01557 

0 

-31070 

0 

0 

0 

0 

-.01557 

P(l.Ole) 

18028 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

260078 

0 

l 

0 

0 

0 

0 

0 

0 

1 

1 

1 

u(L.Olc) 

2.24 

2.1  x  106 

1.05  x  10 

. 

1.55  x  105 

5l 

3.25  x  10 

-206 

0  . 

0 

0 

-.01557 

0 

-3244 

0 

0 

0 

0 

-.01557 

P(l. 02e) 

18679 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

270062 

0 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

u(i.02g 

2.18 

2.12  x  106 

1.06  x  106 

1.06  x  105 

3.29  x  105 

8.17  x  105 

TABLE  14 


EXTRAPOLATED  SOLUTIONS  VERSUS  CALCULATED 
SOLUTIONS1  FOR  PROBLEM  P(e) 


P(0.95e) 

P(0.98e) 

P(0. 99c) 

P<£) 

P(1.01e) 

P(l. 02e) 

P(1.05e) 

55.39 

56.84 

- 

58.30 

-  58J9  _ 

60.25 

63.29 

61.54 

57.82 

58.39 

58.97 

61.52 

_5l-33_ 

56.82 

_ . 

58.31 

.  58^81  _ 

60.30 

63.29 

61.54 

i 

57.82 

58.39 

58.97 

61.52 

Ig99| 

mgm 

wmm 

- 

.  36^23 

36.78 

38.42 

WBEm 

1  s 

35.69 

36.19 

38.39 

1.0 

1.03 

1.04 

__  _ 

1.06 

1.07 

1.10 

1.0 

1.03 

1.04 

1.05 

1.06 

1.07 

1.10 

_  1.0  _ 

1.03 

1.04 

- 

.  __!^06  _ 

JL07  _ 

1.10 

1.0 

1.03 

1.04 

1.05 

1.06 

1.07 

1.10 

1.0 

1.03 

1.04 

__  _  _ 

__1^06  _ 

_  JL07.  _ 

1.10 

1.13 

1.12 

1.04 

1.05 

1.06 

1.07 

1.11 

Entries  have  been  rounded  to  the  nearest  two  decimal  points. 

2 

e,  extrapolated 
3 

c,  calculated 
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As  suggested  by  Figure  LI,  tlie  extrapolated  results  for  Problems 
P(0.99i)  and  P(1.02r)  are  virtually  the  same  as  those  obtained  by  direct 
calculation,  while  the  results  for  Problems  P(0.95t  )  ,  P(0.98c)  and 
l'(l.05i)  arc>  quite  different.  For  example,  linear  extrapolation  from  P(s  ) 
to  P(0.9r>t  )  underestimates  variable  by  about  6  cm  ,  equivalent  to  about 

19.3  percent  of  its  correct  value.  (See  Table  14). 

It  was  also  attempted  to  estimate  the  penalty  function  solution 
"trajectory"  x(>,r)  for  the  perturbed  Problem  P(f:),  using  the  solution  point 
and  corresponding  sensitivity  information  along  the  solution  trajectory  of 
Problem  P(t).  As  shown  in  Table  15  for  subproblem  2,  with  penalty  function 
parameter  r=100  ,  the  results  of  this  extrapolation  for  e  between  .  99> 
and  1 . 02t:  are  quite  satisfactory. 

Figure  12  indicates  the  stability  of  f*(t),  Xj(t),  u^r)  and  their 

partial  derivatives  with  respect  to  the  parameter  i  ^  along  the  solution 
trajectory  of  Problem  P(1.05r).  As  seen,  the  values  for  these  quantities  are. 
again  fairly  stabilized  from  the  third  subproblem  onwards. 
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Figure  11.  x^(e),  XjCe),  x^(e)  and  x^(e)  versus  8 

for  Problem  P(e),  where  e  ■  ^  (extrapolated 
compared  to  calculated). 


3 


TABLE  15 

EXTRAPOLATED  SOLUTIONS  x(c,r)  VERSUS  CALCULATED 
SOLUTIONS1  FOR  SUBPROBLEM  21  OF  PROBLEM  P(f  ) 


P(0.95£) 

P(0.98£) 

P(0.99£) 

P<£) 

P(1.0l£) 

P(1.02f ) 

P(1.05e.) 

H 

s 

55.37 

56.84 

57.34 

58.32 

58.81 

60.29 

1 

1 

63.02 

61.53 

57.25 

57.82 

- - 

58.42 

59.01 

61.64 

e 

55.33 

56.83 

57.34 

58.34 

58.85 

60.36 

2 

c 

63.04 

61.53 

57.26 

57.82 

58.43 

59.02 

61.65 

e  _ 

32.95 

34.60 

35.14 

36.24 

36.79 

38.43 

3 

c 

38.93 

35.90 

35.19 

35.69 

36.20 

36.71 

i 

1 

_  1.0 

1.02 

1.03 

- 

_ 1.08  _ 

1.12 

B 

■ 

1.0 

1.03 

1.04 

1.05 

1.07 

1.10 

xr 

e 

.99 

1.04 

_  _ 

1.06 

1.07 

1.11 

5 

c 

1.0 

1.03 

1.04 

1.05 

1.06 

1.07 

1.10 

X, 

e  _ 

_  J.. 0  _  _ 

1.03 

1.04 

__  _ 

_  I- 06_ 

1.07 

1.10 

6 

c 

1.13 

1.04 

1.05 

1.06 

1.07 

1.12 

All  entries  are  rounded  to  the  nearest  two  decimal  points. 

2 

The  penalty  function  parameter  r  =  100  in  subproblem  2. 

3 

e,  extrapolated 

4 

c,  calculated 
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(k  ■  Subproblea  busbar) 

Figure  12.  Staple  of  stibproblea  solution  point  and  partial 

derivatives  along  the  trajectory  for  Problem  P (1.05c). 
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6.  FURTHER  OBSERVATIONS  AND  CONCLUSIONS 

As  was  noted,  one  of  the  principal  uses  of  sensitivity  analysis  is 
the  estimation  of  a  soLution  of  a  perturbed  problem,  based  on  a  solution  and 
sensitivity  information  of  the  unperturbed  problem.  Obviously  the  "validity" 
of  the  extrapolated  results  depends  on  the  magnitude  of  the  perturbations 
introduced  and,  as  true  even  in  linear  programming,  extrapolations  are  gen¬ 
erally  valid  only  for  a  "limited"  range  of  perturbations.  A  crucial  differ¬ 
ence  is  the  fact  that  sharp  bounds  on  this  range  can  usually  be  determined  in 
linear  programming,  while  this  is  not  usually  true  in  nonlinear  programming. 

In  fact  the  state  of  the  art  is  such  that  we  are  still  developing  techniques 
for  calculating  sharp  error  bounds  for  a  specified  measure  of  stability,  an 
objective  being  to  calculate  such  bounds  in  a  computationally  efficient 
manner.  It  is  not  satisfactory  to  require  as  much  effort  to  calculate  a  sensi 
tivity  measure  as  it  takes  to  solve  the  given  problem  from  scratch.  This  is 
often  not  feasible  in  practice,  nor  does  it  appear  to  be  generally  essential 
in  theory.  A  whole  body  of  methodology  has  yet  to  be  developed  before  useful 
stability  calculations  "in  the  large"  are  rigorously  validated  and  efficiently 
implemented  for  general  parametric  nonlinear  programming  perturbation  analysis 

Based  on  the  results  of  the  present  analysis,  the  following  conclusions 
are  also  offered,  subject  to  the  scenario  described  herein: 

1.  The  weight  of  the  corrugated  bulkhead  is  extremely  sensitive 
to  changes  in  the  minimum  allowable  panel  thickness,  which  is 
in  turn  a  function  of  the  ship  length.  (This  latter  relation¬ 
ship  is  not  present  in  the  given  model  or  pursued  in  this  paper.) 

One  millimeter  relaxation  of  this  parameter  will  save  about  0.331 
tons  of  steel  per  bulkhead  (equivalent  to  6.3  percent  of  optimal 
bulkhead  weight) . 

A  similar  result  holds  for  the  corrosion  allowance  parameter  k^  . 
One  millimeter  relaxation  of  this  parameter  would  save  about 
0.11  tons  of  steel  per  bulkhead  (equivalent  to  about  2.1  percent 
of  the  optimal  bulkhead  weight). 
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2.  Weight  of  the  corrugated  bulkhead  is  only  partially  sensitive 
to  changes  in  the  top  and  middle  panel  lengths  and  total  panel 
depth.  These  sensitivities  are  completely  overshadowed  by 
those  for  minimum  allowable  thickness  of  the  panels. 

3.  Surprisingly,  the  weight  of  the  bulkhead  is  rather  insensitive 
to  changes  in  the  parameter  k^  (a  function  of  the  allowable 

2 

bending  stress).  A  50  kg/cm  relaxation  on  maximum  allowable 
bending  stress,  which  is  a  rather  significant  relaxation  as  far 
as  design  considerations  are  concerned,  would  save  only  about 
0.0162  tons  of  steel  per  bulkhead  (about  0.3  percent  of  optimal 
bulkhead  weight) . 

4.  Assuming  the  validity  of  the  model,  these  results  suggest  that 
the  most  economical  marginal  reductions  in  bulkhead  weight  will 

be  realized  by  first  decreasing  the  minimum  plate  thickness  require¬ 
ment  (tmin)  and  next  by  relaxing  the  corrosion  allowance  (k2)  . 

5.  For  rather  significant  perturbations,  the  extrapolated  solution 
of  the  perturbed  problem,  based  on  the  solution  and  sensitivity 
information  of  the  unperturbed  problem  is  very  accurate  for  the 
subject  problem. 

As  a  concluding  remark,  it  is  universally  acknowledged  that  a  sensitivity 
analysis  is  an  indispensable  part  of  a  "solution,"  certainly  in  any  practical 
context.  Aside  from  providing  invaluable  insight  into  problem  and  model 
structure  and  estimates  of  the  effects  of  changes  in  design  or  data  parameters, 
a  thorough  analysis  can  provide  warnings  of  instabilities,  indications  of 
improper  problem  or  model  formulations,  and  detailed  guidelines  for  making 
cost-effective  changes  in  the  parameters. 

The  computer  program  SENSUMT  is  available  to  the  general  public  on 
request  to  the  authors. 
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